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Unconcatenated, unknotted polymer rings in the melt are subject to 
strong interactions with neighboring chains due to the presence of topological 
constraints. We study this by computer simulation using the bond-fluctuation 
algorithm for chains with up to = 512 statistical segments at a volume 
fraction $ = 0.5 and show that rings in the melt are more compact than 
gaussian chains. A careful finite size analysis of the average ring size R oc N'^ 
yields an exponent v ~ 0.39 it 0.03 in agreement with a Flory-like argument 
for the topological interactions. We show (using the same algorithm) that 
the dynamics of molten rings is similar to that of linear chains of the same 
mass, confirming recent experimental findings. The diffusion constant varies 
effectively as oc A'^~i-^^(^) and is slightly higher than that of corresponding 
linear chains. For the ring sizes considered (up to 256 statistical segments) 
we find only one characteristic time scale Tee oc N'^-'^C^^-^ this is shown by the 
collapse of several mean-square displacements and correlation functions onto 
corresponding master curves. Because of the shrunken state of the chain, this 
scaling is not compatible with simple Rouse motion. It applies for all sizes of 
ring studied and no sign of a crossover to any entangled regime is found. 

PACS numbers: 61.25.Hq, 61.41. -Fe 



I. INTRODUCTION 

Ring polymers have been extensively studied experimentally by several groups, not al- 
ways leading to completely consistent results l^Hlfljll • Indeed, the synthesis of the samples 
of unknotted and nonconcatenated rings (see fig. |I|) remains a very delicate issue 0. Nev- 
ertheless, the general conclusion emerging from these studies is that the dynamics of rings 
is quite similar to that of linear chains of the same molecular mass and density . This is 
certainly a surprising result from the point of view of the reptation model which describes 
the flow behavior of entangled linear polymer chains @,T^]. In this effectively single-chain 



model, the constraints on the motion of a reference chain due to the interactions with its 
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neighbors are replaced by a curvilinear "tube" within which the chain "reptates" : relaxation 
of the constraints occurs only at the chain ends. Clearly, a closed ring polymer (having no 
ends) cannot reptate in this conventional sense. Accordingly, the motion in a melt of rings 
should be quite different from that of an analogous linear chain system; and at first sight 



on would expect it to be slower. Indeed it was argued [|ri|,|Ty] that motion of ring polymers 
should be exponentially slow and comparable with that of star polymers fl^ . 

This is however not borne out by rheological measurements on polystyrene (PS) and 
polybutadyene (PB) rings, which showed that the zero-shear viscosities tjq for melts of rings 
for all molecular weights considered are similar to, but even slightly smaller than, those for 
linear chains The temperature dependence of 770 in PS rings is virtually indistinguish- 
able from that of linear polystyrene of high enough molecular weight; all viscosities display 
classical WLF or Vogel type temperature dependence and they can be superimposed on a 
single master curve. For linear chains, two power-law regimes arise in the mass dependence 
of the viscosity, corresponding to Rouselike and entangled motion. Unexpectedly, the same 
is found for ring polymers. Below a critical mass the viscosity seems to increase with an 
exponent ai ^3/2 which is larger than the expected ai = 1 for ideal Rouse behavior. The 
values for rings and linear chains of the same molecular weight are of similar magnitude. In 
the high mass regime (M > M^) the viscosities increase more strongly, again with a power 
similar to the linear case. An exponential increase in the range of masses available (up to 
M ^ 1.85 ■ 10^) is explicitly excluded gig. 




FIG. 1. Sketch of the opposed topological constraints, a) On the left is a permitted configuration 
of an unknotted ring. It cannot turn into the forbidden (knotted) configuration on the right (nor 
vice versa), b) On the left is a permitted configuration of a pair of unknotted, unconcatenated 
rings. It cannot turn into the forbidden (concatenated) configuration on the right (nor vice versa). 

The similarity of the dynamics of ring molecules to their linear counterparts is also re- 
ported on a more microscopic level from tracer diffusion measurements 0,^]. Dilute labeled 
PS chains of mass M, in a matrix of mass P, were measured using Forward Recoil Spec- 
troscopy. Different topologies (rings in linear matrices [0] , linear chains in ring matrices [§ , 
rings in microgel [^) were investigated and compared to linear chains in linear matrices. 
Unfortunately, there exist so far no systematic measurements of ring tracers in matrices of 
rings (of identical size) in which the molecular mass covers a significant range. However, 
consistent with the rheological measurements quoted above, the tracer diffusion of linear PS 
in ring matrices was found to be nearly identical to that of linear PS in linear PS matrices 



T5[| . This is again a surprising result in view of theoretical concepts of matrix dependent 
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tracer diffusion |TT[ . 

A suggestion by Lodge et al. |10| is that the experiments on rings may be reconciled with 
the reptation concept by talcing into account the higher entanglement length. The critical 
mass for entanglement can be estimated from the viscosities or from the shear modulus. 
PS rings exhibit a plateau modulus approximately one-half (one-fifth for PB) that for linear 
polymers, suggesting that rings are less effective at forming entanglements. Rubber elasticity 
arguments indicate a critical mass for polystyrene cycles of Mc = 58000 compared to Mc = 
30000 for linear chains. If one now makes a comparison at an equal number of entanglements 
ring melts in fact have higher viscosities than linear chains. The experimental viscosity data 
then lie in the unentangled/entangled crossover region, where it is clear from the linear 
polymer results that reptation is not the only available mode. This explanation would allow 
a crossover to much slower ring dynamics at very high molecular weights. 

No theoretical explanation for this increase of the critical mass exists, but it seems 
reasonable to relate it to the idea that ring conformations in the melt may be partially 
collapsed Besides the usual excluded volume interaction between two neighboring 



chains there is a topological interaction due to the exclusion of knotted and concatenated ring 
configurations (see fig. |l|). While the excluded volume interaction is screened out whenever 
the overlap with other chains is large (i.e., as long W^/N » 1, true for large chains if 
u > 1/d) the topological interaction in 3-dimensions need not be screened and should cause 
a reduction in the size exponent u below the gaussian value (z^ = 1/2). 

Despite this argument, no experimental study of the radius of gyration of rings in the 
melt appears to have been made. A good understanding of the static properties of the system 
is of course an indispensable starting point for a reasonable description of the dynamics, so 
this is unfortunate. Detailed static measurements (in both dilute and concentrated systems) 
could also provide a stringent test on the quality of the ring synthesis, as McKenna et al. 
have pointed out 0. Without them, for many synthesis routes it is hard to exclude the 
possibility of knotted or concatenated rings, nor that of a small fraction of linear chain 
contaminant which might modify the viscosity substantially (by threading the rings, for 
example) 

In view of those experimental difficulties, a detailed simulation study to investigate both 
the statics and the dynamics in a melt of rings is certainly warranted. Most previous at- 
tempts are based on Pakula's "cooperative motion algorithm" p!7| , p!8[| . While this algorithm 



can give correct results for the statics (especially for melts) its value for studying dynam- 
ics is questionable (see section 2 below). In the following we study melts of non-knotted, 
non-concatenated rings within the well-established "bond-fluctuation model" (BFM) to get 
a better insight in their statics and dynamics. Apart from its computational efficiency, an 
advantage of the method is that comprehensive data for linear chains have already been 
obtained with which quantitative comparisons can be made. Indeed, even if ring 

polymers did not exist experimentally, this comparison might help illuminate several long- 



standing but still controversial issues in the dynamics of entangled linear chains |21 



The article is organized as follows : In the next section we give briefly some technical 
comments on the BFM simulation performed. In section 3 results on the conformation 
of isolated rings (in an athermal solvent) are presented, showing that the unknottedness 
constraint is insufficient to swell the rings significantly beyond what would be caused by 
excluded volume forces alone. Molten rings, on the other hand, are found to be quite 
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compact with an exponent v 0.4, consistent with the assumption [|T6| that roughly one 
degree of freedom is lost for every topological interaction with a neighboring chain. We then 
show, in section 4, that the dynamics of rings and linear chains are qualitatively similar 
over much of the range of chain length that we are able to simulate. For the largest rings, 
however, a perfect scaling (involving a single characteristic time for each chain length) is 
still obeyed whereas for this size of linear chain significant departures are observed. This 
allows us to confirm the suggestion that the entanglement mass is much larger for rings. 

We also study the degree to which rings thread through one another, a question addressed 
in both sections 3 and 4. This seems to be insufficient to allow large clusters of mutually 
entangled material to build up. The extent of entanglement, quantified roughly as the 
number of neighboring rings contacting a given molecule, nevertheless appears relevant for 
dynamics: for both rings and their linear counterparts the diffusion constant is shown to 
scale as the mass of the correlation hole (at least, in the range of chain lengths studied). 
Finally in section 5 we give our conclusions and discuss the impact of these ring investigations 
on our understanding of entangled polymer dynamics. 



II. A BOND-FLUCTUATION MODEL STUDY 

The algorithm used in this investigation is the well-established "bond-fluctuation model" 



(BFM) of K. Kremer et al [g^. This coarse-grained 3-dimensional lattice model has proven 



to be especially useful for investigating the universal features of statics and dynamics in dense 



polymeric melts |jT9[. A small number of chemical repeat units (i.e. a Kuhnian segment) 
is mapped onto a lattice monomer such that the relevant characteristics of polymers are 
retained: connectivity of the monomers along a chain and excluded volume of the monomers. 
Each monomer occupies a whole unit cell of a simple cubic lattice with periodic boundary 
conditions. Adjacent monomers along a polymer are connected via one of 108 allowed bond 
vectors. These are chosen such that the local excluded volume interactions prevents the 
chains from crossing each other during their motion. This conservation of the topology 
ensures that the rings (which are set up initially as thin "sausages" ) remain neither knotted 
with themselves nor concatenated with one another during the relaxation and sampling. In 
our athermal simulation, "local jumps" are realized by choosing one monomer at random 
and attempting to jump over the distance of one lattice spacing in one of the six basic 
directions (also randomly chosen). The attempt is accepted if excluded volume restriction 
are satisfled and the new bond vectors to the neighbors along the ring belong to the allowed 
set. 

There has already been made a simulation study of ring polymers by Frisch et al in 
the framework of the BFM. However, due to the large CPU-time demands, that study was 
restricted to single rings; also the dynamical properties were not very directly addressed. 
The only existing simulation data on the statics and dynamics of rings in the melt use 
Pakula's "Cooperative Rearrangement Algorithm" [p4| , ^ . In contrast to the local jumps 
utilized in the BFM, this algorithm changes simultaneously and collectively the monomer 
positions on a number of different chains. While giving correctly the static properties, a 
clear correspondence between Monte-Carlo time and real time is yet to be established for 
this algorithm and its dynamical interpretation is accordingly unclear. 

In the present BFM investigation we want to extend the careful study of linear chains 
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made by Paul et al to ring polymers and compare our results to their data on linear 
chains. At a filling fraction $ = 0.5 of occupied lattice sites, many static and dynamic 
features of molten polymeric materials are reproduced by the BFM. For example, the single 
chain conformations obey gaussian statistics down to the screening length ^ ^ 6 (in units of 
the lattice constant) of the excluded volume interaction obtained from the static structure 
factor [l^ . There are extensive results on the dynamical properties covering the range from 
an unentangled behavior for short chain lengths up to the onset of reptation-like motion for 
chain length = 200. Of course, the dynamics of long polymers in a dense melt slows down 
dramatically with growing chain length and therefore poses huge demands on CPU-time re- 
quirements. For the present investigation we employ a very efficient implementation of the 
BFM on a massively parallel CRAY T3D [^. Using a 2 dimensional geometrical decom- 
position of the simulation grid of linear extension L = 128, we employ 64 T3D processors. 
This permits us to equilibrate systems comprising 131072 monomers and ring lengths up to 
512 (statics) or 256 (dynamics) statistical segments. This study involved about 5000 hours 
of single processor CPU time. 

The starting configurations consisted of straight sausage-like rings (loops enclosing no 
area) which were carefully equilibrated for at least one relaxation time (i.e. the center 
of mass had moved a distance comparable to the chain size) before any data were taken. 
Indeed, for all but the longest chains, runs were continued well beyond this (so as to generate 
dynamical data) which enabled us to confirm that the static chain extensions had settled to 
their equilibrium values by this time. 



III. STATICS : CONFORMATIONS OF RINGS IN THE MELT 

In this section we consider the effects of the unknottedness and nonconcatenation con- 
straints on the conformational properties of dilute and molten rings. While the effect of 
the former constraint on isolated rings turns out to be irrelevant the non-concatenation 
requirement significantly compactifies molten rings. 

The size of the rings is measured firstly with the usual mean-square radius of gyration 
Rfj] as a second measure, we define the average distance between pairs of monomers 

that are A^/2 monomers apart along the ring contour, {Rl) = (^(^Rn — Rn+N/2^ ^, and call 

this the mean-square ring diameter. For isolated rings in an athermal solvent results for 
both quantities are presented in table |. In agreement with previous Monte Carlo studies 
by Frisch et al |2^ the simulation yields a ring size R oc with exponent u ~ 0.595 



for the radius of gyration and u ~ 0.605 for the ring diameter. These values are only 
slightly larger than the excluded volume exponent (z/ = 0.588) for linear athermal chains, 
the difference lying within the range of the statistical error. The influence of the topological 
constraints on the static properties of isolated rings thus appears not to alter the chain- 
swelling exponent in 3 dimensions (though there may be a prefactor effect) ||2^. This 
flnding concurs with the analytical studies of des Cloizeaux [^] and is now also conflrmed 
experimentally. Indeed, while early SANS data of J. Higgins et al [H seemed to indicate 
a nearly gaussian statistics for rings in good solvent, later the dehcate dependency of the 
ring properties on preparation conditions was overcome by G. Hadziioannou et al [Q] and 
J. Roovers et al corroborating that the statistics are as for linear chains in a good 
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solvent. However, a decrease of the ^-temperature of isolated rings (compared to linear 
chains) by several Kelvin has been measured. This is defined by the point at which the 
second virial coefficient from the light scattering measurements equals zero, and seems to be 
the only manifestation of the unknottedness constraint A Monte Carlo study of the 

0-point depression of isolated rings, including a careful finite-size analysis could certainly 
yield interesting additional information; we do not attempt this here. 

The situation changes completely in the other limit of molten rings where, confirming 
Pakula's simulation we find very compact rings (table H). When naively fitted without 
any finite size analysis, the ring sizes yield an exponent v ^ 0.44 for the ring diameter 
Re and a slightly higher value, u ~ 0.45, for the radius of gyration Rg. Those values are 
quite similar to the ones obtained by Pakula. However, finite-size corrections are clearly 
detactable in the curvature of the data points, and a more careful fitting is necessary, as 
explained below and in fig. This yields in an extrapolation to the limit of infinite masses 
a distinctly lower exponent u = 0.39 ± 0.03, which (as it should) becomes the same for our 
two measures of ring size, R^ and Rg. 
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FIG. 2. The running exponents i^(A^) obtained from the radius of gyration \^g^ (circles) and 
the ring diameter (-Rg) (squares) versus the natural variable of the screening . Both the two 
exponents and the ratio (iie) / (^g^ the inset lay on straight lines. From this finite size analysis 
we obtain the exponent v ~ 0.39 it 0.03 in the limit of infinite mass. 

At the volume fraction $ = 0.5 used in the simulation, the excluded volume interaction 
is not completely screened, as pointed out in the previous section. This explains why the 
radius of gyration, which is more sensitive to short-scale structure, tends to give larger 
exponents z^(A^) than the diameter which probes larger distances (fig. |^). At high masses 
this short-scale effect of the excluded volume will become less and less important, so both 
the ratio of diameter to the radius of gyration (inset of fig. |^), and the running exponents 
z/(A^), tend to become smaller. These exponents are plotted against the natural variable 
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of the screening (proportional to i/Rg) [^H; the procedure brings all the measured values 
z/(A^) to lie on straight lines which permit a precise finite size scaling analysis [^|. As 



one sees from fig. ^ the low masses gives exponents close to those obtained previously; the 
highest masses used yield an exponent v ^ 0.42. Only in the limit of infinite masses does 
the short-scale excluded volume effect vanish and the two lines for Rg and Rg merge at a 
value u = 0.39 ± 0.03. 

This value is consistent with a crude but interesting Flory-like estimate of the free energy 
of a polymer ring given by Gates and Deutsch [0]. They argued that, if a ring of polymer- 
ization index N has size R, it is overlapped with a number of neighboring rings of order 
R'^/N (in d dimensions). The more spatially extended the ring, the more entropy is lost by 
the nonconcatenation constraint with its neighbors. The simplest possible estimate of this 
is to say that the number of degrees of freedom lost due to the constraint is proportional 
to the number R'^/N of neighbors which the ring is prevented from threading. This gives 
a contribution to the free energy of F oc kTR'^/N tending to decrease the ring size. On 
the other hand, there is also an entropy penalty if the ring becomes too squashed; the free 
energy required to squash a gaussian chain of steps into a region of linear size R less than 
N^^"^ scales as kTN/R^. Adding these contributions, and minimizing over R, gives a charac- 
teristic size scaling as i? oc N'^ with u = 2/{d + 2). Hence in three dimensions the exponent 
is z/ = 2/5, very close to the value found by our simulation. Note that the latter value is 
definitely smaller than the result u = 1/2 — l/Gn put forward recently by M. Brereton and 
T. Vilgis m. 
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FIG. 3. Distribution function of the ring diameter Rg as function of the reduced variable 
Re/{Rlf'^. Data for various masses are distinguished by different symbols, as indicated in the 
figure. Apart from the very small chains all distributions superimpose. 

We note that from the measured exponents of the chain size of an isolated chain z/q ~ 0.6 
and of ring chains z/^ ~ 0.4 a crossover scaling can be obtained in the usual way. This 
yields for the ring size in the semidilute concentration range, R ^ Ro ($/$*)"^/^ where 
the monomer density is $ and the usual crossover density (overlap threshold) is $*. The 
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decrease in the size of a ring with $ is much more pronounced than for hnear chains where 
R ~ $-1/8. 

In fig. ^, we show simulation data for the probability distribution P{Re) of the diameter 
of a ring as a function of the characteristic variable x = Re/{Rl)^^'^- Data is given for 
various masses up to = 256. The scaling collapse becomes more and more perfect with 
increasing mass indicating again that the local interactions become irrelevant for ring sizes 
much larger than the excluded volume screening length ^. In principle, this distribution 
should define some further characteristic exponents, assuming that, as for linear chains 
3^1 the distribution rises with a power P{x) oc for small x and drops off essentially as 



P{x) oc 2;"exp(— x*^) for large x. For isolated linear chains the exponents g and 6 can be 
written in terms of u and 7 (the latter is the exponent controlling the dependence of the 
free energy). Unfortunately, although the scaling in fig. |^ is good, the data is not precise 
enough to extract any corresponding exponents in this case. Accordingly we leave this issue 
for future investigations. 

A quantity of more direct experimental relevance than P{r) is the static structure factor 
S{q). In fig. ^ the values of S{q)/N for masses up to A^ = 512, expressed as function of the 
characteristic variable RgQ, superimpose (apart from the spurious Bragg peak in each case) 
onto a single master curve. Consistently with our earlier discussion of the scaling of the ring 
size, a self-similar power-law regime is apparent, and shows a fractal dimension 1/z/ = 1/0.4 



as expected. The Gaussian dimension 2 is clearly ruled out by our data |35 




FIG. 4. Structure factor S{q)/N versus the characteristic variable Rgq for various masses A^ as 
indicated in the figure. The dashed line confirms the fractal dimension 1/u = 2.5 = 1/0.4. The 
sohd slope 2 for rings with gaussian statistics cannot match the measured structure factor. 

As mentioned in the introduction, and above, the probability of threading of a ring 
polymer by its neighbors is implicated in the mechanism of partial collapse, and may also 
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have dynamical consequences by way of long-lived clusters of entangled material slowing 
down the relaxation times. It is difficult to define precisely what is meant by threading (in 
either static or dynamical terms) but a partial measure of the ease with which a ring can 
be threaded is to measure the area of its projection onto a random direction. This area a 
is defined as a signed quantity (the component in that direction of the vector area of the 
ring) which vanishes for any configuration in which the ring exactly retraces its own steps. 
A measure of the "threadability" is provided by A = it turns out that A cc Rg as is 

shown in fig. |^. The scaling is the same as would be naively expected. Note that sausage-like 
configurations would involve very small values of A. The magnitude A/R"^ ^ 1 indicates 
that the rings have no tendency to retrace their own steps; a fact also revealed by inspection 
of snapshots. 

Another quantity that may be relevant in measuring the strength of the topological 
interactions is the number of chains touching a reference chain n^y. Two chains are defined 
to be "touching" whenever their monomers include pairs separated by a distance less or 
equal to \/6 lattice constants. (This somewhat arbitrary microscopic distance was chosen 
to include all the monomers within the first peak of the monomer-monomer correlation 
function.) Intuitively njy should vary as the number of chains within the "correlation hole" 
spanned by a given chain, and we see in fig. |^ that un indeed scales as R^/N. The plateau 
at high enough masses confirms the self-similarity of rings in the melt which was already 
indicated by the scale invariance of the diameter distribution in fig. |[ 
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FIG. 5. Surface circumscribed by a ring in the melt A, plotted as A/R'^ (circles) and the number 



of neighbors n^, plotted as nN^/Rg (squares) versus A^. The mass-independent plateaus 0. 
3 respectively are reached for masses larger than = 100. 
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IV. DYNAMICS OF RINGS IN THE MELT 



As pointed out in the introduction, there is some experimental evidence that the dy- 
namics of rings are similar to their linear counterparts, at least up to the largest molecular 
masses that can readily be obtained. Our simulation data, presented below, confirms this for 
rings of up to = 256 monomers by comparing mean-square displacements, chain and co- 
operative motion correlation functions, and the resulting diffusion coefficients and relaxation 
times. 




FIG. 6. Mean-square displacements of a ring monomer gi{t), of a monomer in the frame of the 
center of mass g2{t) and of the center of mass gsit) versus the reduced time t/ree- The three times 
Ti^2,3 (vertical dashed Unes) are defined in the text. The two horizontal dashed Unes correspond 
to mean-square displacements of 2/3 ^ii^^ and (^-^g^ respectively. The Fickian behavior at long 
times is indicated by the solid line. The effective exponent x = 0.81 of the center of mass motion at 
shorter times due to net forces of neighboring chains is displayed by the broken line. The anomalous 
diffusion of a monomer for molten rings with exponent v = 0.4 is indicated by the dashed-dotted 
line. 

We characterize the dynamics by measuring three different mean-square displacement 
functions, describing the motion of monomers gi{t) = ((R„(t) — R„(0))^) in the lab frame, 
the motion of monomers in the center-of-mass frame of a given ring, g2{t) = ((R„(t) — 
R-cml^) — R-n(0) + Rcm(0))^)/2, and the motiou of the center of mass itself, g^it) = {(Remit) — 
R'cm(O))^). The mean-square displacement functions for molten rings are given in fig. |^. 
These quantities are scaled by the mean square characteristic ring size (Rg) whereas the 
time coordinate is scaled by the rotational relaxation time Tee- This rotational time is the 
decay time for relaxation of a ring diameter (the vectorial displacement between monomers 
N/2 apart in the sequence) which is obtained from the correlation function Cee{t) described 
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below (shown in fig. 0), and its scaling with is shown in fig. ([TT|). Apart from the mass 
= 16 all curves perfectly superimpose with the use of this single scale factor for each chain. 
This indicates that no second time-scale is present as would be expected for an entangled 



system (the entanglement time). Such scaling is sometimes called Rouse behaviour [jT9 
but here we describe it as "unentangled scaling" since the power laws involved need not be 
those of the Rouse model. Notably, no crossover whatever is detectable from this scaling to 
a modified (reptation or other) motion, involving a second time scale, even for the largest 
rings studied (A^ = 256). This contrasts with the fact that for linear chains of mass larger 
than A^ = 200 the unentangled scaling starts to break down, with clear signs of a crossover 
to a new regime (whose exponents could not, however, be reliably measured) This 
finding confirms that any entanglement length for rings is larger than that for linear chains 
of the same type, in agreement with experimental inference 

As mentioned before, it is tempting to associate this with the fact that the rings are 
partially collapsed. Crudely one could argue that for a given degree of entanglement, the 
number of chains in the correlation hole (which scales as Rg/N) should be the same for 
rings and linear chains. This would give a rough estimate of A^ = 1000 BFM monomers as 
the breakdown point of the unentangled scaling in the ring case. This estimate is broadly 
consistent with the entanglement mass for rings reported in experiments as lying between 2 
(for PS) or 5 (for PB) times as large as for linear chains, as discussed in the introduction. 

Qualitatively the mean square displacements of rings displayed here show very similar 



behavior to that of short linear chains, up to masses of order A^ = 100 [p!9| , p0[ . (However as 
mentioned above, for larger masses the scaling breaks down in the linear chain case, and in 
this regime the rings and linear chains are no longer precisely alike.) While the center of mass 
follows the Fickian type of diffusion at long times (t/ree 3> 1) a clear signature of net forces 
acting on the center of mass of the chains is displayed at short times, when the center of mass 
displacement is proportional to (t/ree)^, with an effective exponent x ~ 0.81. (Within a 
true Rouse model, in which each monomer in the system is subject to uncorrelated random 
forces, this exponent would have the Fickian value of unity at all times.) This effect was 
also observed for linear chains with a slightly higher effective exponent x ~ 0.85 at the same 



density $ = 0.5 ||T9|,^. There it was also verified that the exponent x is density dependent, 
approaching 1 in the low density limit; it is very likely that this is similar for rings. We 
can therefore conclude that, as for linear chains, there is a net force acting on the center 
of mass generated by the interaction of the test chain with surrounding ones, slowing down 
the chain motion at short times. 

The monomer motion may also be split into two regimes characteristic of short and 
long times. While in the latter the monomeric displacements approach asymptotically that 
of the whole chain {gi{t) oc t), the conformational properties of the rings in the melt are 
refiected in the short time anomalous diffusion regime. The involvement of more and more 
ring monomers as the lifetime of a fluctuation increases, gives by general scaling arguments 
a mean-square monomer displacement of gi oc (t/reeY^^^^^^'^''^ ■ For our rings (z/ «i 0.4) 
this yields an exponent of about 0.45 which agrees well with the dynamical simulation data 
shown in fig. ^. (The agreement is even better if one uses for each chain length the measured 
exponent i^(A^) which is somewhat larger than the extrapolated value z^(oo), as explained in 
the previous section.) Defining a monomeric mobility W hj gi = b'^{WtY^^-^^^^'^'^^ where b is 
the monomer size, perhaps surprisingly we find (as shown in tab. |I| ) that W ^ 7 ■ 10~^ for 
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rings is about 4 times larger than for linear chains where W ^ 1.6 ■ 10~^. The mean-square 
monomer displacement in the frame of the center of mass, g2, follows (as expected) for small 
times the same behavior of that in the lab frame gi while for long times it approaches (by 
definition) the mean-square radius of gyration of the rings. 

In order to characterize succinctly the mean square displacement curves for rings, and to 
compare them to the linear chain counterparts, we define (following ref. [jl9l) three charac- 
teristic times, ri_2,3 according to the following criteria: 5'i(ti) = (^Rf), 92(^2) = 2/3 (-R^) and 
g2{jz) = Qzij^)- Each r-parameter is a measure of the decay of the corresponding displace- 
ment function, as indicated in fig. ^ Apart from a bit of scatter we obtain, independently of 
mass, the ratios T2/T1 3.3 (compared to ~ 0.9 for linear chains) and t^/ti 8 (compared 
to ~ 3). This means that, normalizing by the motion of a monomer in the lab frame, both 
the center of mass motion, and that of monomers in the center of mass frame, take longer 
to reach their long time asymptotic limits than is the case for linear chains. 
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FIG. 7. Distribution P[x) of the mean-square displacements x of the center of mass 53 (t) for 
chains of mass = 256. Data over a large range of time (as indicated in the figure) collapse 
perfectly on one straight line with slope —3/2 when plotted as g^^"^ {t)P{x) / x^/"^ versus x/g^{t), 
corresponding to a gaussian distribution of the displacement variable. 

We now discuss the probability distribution P{x) for the mean square center of mass 
displacement x = (Remit) — Rcm(O))^. In linear chain systems, ah chains behave roughly 
alike and, on the time scale of motion over one or more gyration radii, the distribution of 
Rem is essentially gaussian. This would lead to P{x) oc gs{t)~^^'^x~^^'^ exp[—3x/2gs{t)]. For 
rings, another scenario is possible, in which at a given time a small number of rings are 
relatively "unentangled" (for example with wormlike configurations), as discussed by Klein 
|1 1| , while others form entangled clusters which can scarcely move. This would yield, in 
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the crudest picture, a bimodal distribution for P{x) at times shorter than the hfetime of a 
cluster. This possibihty is apparently ruled out by the distributions P{x) we obtained as 
for instance in fig. |^ for a melt of rings of mass N = 256. The distribution for times ranging 
from much smaller than the relaxation time T^e = 1-6 ■ lO^MCS, up to times much larger, 
could be superimposed on one single master curve, which shows precisely the form expected 
for gaussian scaling as considered above. 

From the mean-square displacement of the center of mass gsit) one obtains the diffusion 
coefficients as shown in fig. ^ for isolated rings and in fig. ^ for rings at a volume 
fraction $ = 0.5. For isolated rings in athermal conditions, the diffusion constants obey 
D]\r ^ 0.034/A^; these are virtually identical to their linear counterparts (fig. |]), as one 
expects from an essentially Rouse dynamics in this limit (hydrodynamic forces are of course 
excluded in our model). The situation is somewhat different for the diffusion constants Dn 
obtained for rings in the melt. The first striking point is that the rings diffuse faster than 
linear chains of same mass and density. This confirms the trend found experimentally in 
zero-shear viscosities for PS and PB melts of rings [^]. 
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FIG. 8. Diffusion coefficient ND^ versus mass for isolated rings (circles) and linear chains 
(squares) . 
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FIG. 9. Diffusion constant of rings (circles) and linear polymer (squares) as a function of the 
mass A^. The diffusion constant of rings varies effectively as ND]\f oc N~^''^'^ (solid line). Neglecting 
any curvature of the data points the diffusion constants of linear chains are shown (dashed line) to 
vary as NDj\f oc N^'^. This is compared to the reptation prediction for linear chains (dashed-dotted 
line). 

Second, for rings, there is no sign of curvature on the log plot, which is consistent with 
the idea that the entanglement mass for rings is larger than for linear chains. In fact, the 
diffusion constant of linear chains, also shown in fig. has previously been interpreted in 
terms of a crossover from true Rouse behaviour for small N (which would appear as a plateau 
on this representation) to a new, entangled regime, with the crossover effects first appearing 
for 100 < N < 200. For linear chains (but not rings) this crossover behaviour was much 
more strongly seen in gi and g^. Looking at the center of mass diffusion data alone, however, 
there is no sign of a small- plateau, and no firm evidence of a crossover to a new regime of 
entangled behaviour, even in the linear chain case! This suggests that other interpretations 
might also be worth investigating. 

If one assumes that, within the range of masses studied, both linear chains and rings 
follow a single power law behavior, then one finds respectively D ~ A^~^'^ and D ~ A^~^'^ 
for the two cases |3^. Remarkably therefore, in each case that the diffusion constant varies 
inversely with the mass of the "correlation hole" of surrounding chains: D]\f oc l/Rg- This 
suggests a mental picture in which the center-of-mass mobility of any given chain is governed 
by having to "drag" the contents of the correlation hole along with it. Although this picture 
should clearly break down for long linear chains (where a more conventional reptation pic- 
ture becomes appropriate) it could provide some insight into the behaviour at intermediate 
A^. For rings, this "intermediate" regime appears to be more strongly developed. If one 
compares rings and linear chains at equal mean-square gyration radius, then the rings have 
slightly smaller diffusion constants than linear chains. Assuming that any increase in the 
entanglement length of rings is due to partial collapse, this is consistent with the experimen- 
tal fact that the viscosities of ring melts, compared with linear chains at an equal number 
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of entanglements, have slightly higher viscosities. 

Master curves for some further correlation functions are presented in fig. |10[ The corre- 
lation function Cee(t) is analogous to the end-end vector correlation function usually used 
for linear chains 0. It describes the decay of the diameter vector (Re(t) ■ Re(0)) between 
two monomers of a ring separated by N/2 monomers. The relaxation time Tee, used to scale 
the time axis in the various plots discussed already, is defined as the time at which this 
correlation function has decayed by a factor 1/e. Note that this correlation function shows 
a near-exponential decay. 
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FIG. 10. While the correlation function of the diameter vector Cee{t) (data on the left side of 
the figure) drops off exponentially, the correlation function Cn{t) measuring the mean number of 
chains touching a given reference chain decays much more softly with a power law. Data points 
for various masses superimpose when plotted versus the reduced time t/Tee- 

This contrasts with the second correlation function Cn{t), which measures the decay in 
the mean number of chains, "touching" (in the sense defined above) a given reference chain 
at time zero, that are still touching it (or touching it again) at time t. From this quantity is 
subtracted an asymptotic plateau value (arising from the finite size of the simulation cell) 
so that the correlation function vanishes at long time; we normalize to C„(0) = 1. This 



Cn{t) apparently decays with a power law (fig. [To|). This correlation function Cn{t) was 
defined in an attempt to monitor any possible clustering of entangled rings (leading to a 
fraction of slow- moving material), by detecting possible long-lived contacts between rings 
that might arise from threading of one ring by another. Although the number of contacting 
rings remains large for times much longer than the measured diffusive relaxation times, there 
is nothing like a plateau; the power-law decay of Cn apparently scales with T^e, so there is 
no new time scale due to clustering It seems that the threading probability is not large 
enough to form large clusters whose percolation (for example) could lead to new dynamical 
phenomena. 

In fig. |Tl| we compare the relaxation times Tee as defined from the decay of the diameter 
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correlation function Cee{t) (and used above for the scaling of the time axis) with the quantity 
{R^g) / D which also defines a characteristic time. To a good accuracy, these are proportional, 
but the latter is larger by a factor of about 10. The scaling of the relaxation time, by either 
definition, follows approximately an A^^ law (as do the times ri, T2 and we obtained from 
the mean square displacements), as predicted from the simple Rouse model for gaussian 
chains or rings. However, we believe that this is fortuitous, the diffusion constant being the 
more fundamental quantity; indeed, true Rouse motion for partially collapsed objects (as 
we know the rings to be) would give a higher power. For comparison, the relaxation time 
obtained from the diffusion coefficient via Rg/Dj^ for linear chains is also plotted. In the 
range of studied, this follows a A^^'^ power law. (The linear chain values are taken out of 



Ref. [191, some data for short chains are added.) 



10° 



10' 



10° 



10= 



10^ 



10° 



10 



slope 2.58 




100 
N 



1000 



FIG. 11. Reorientation relaxation time Tee for rings (circles) and linear polymers (squares) 
obtained from Cee{t) versus polymer mass N. The solid line indicate the exponent 2.06 found for 
the collapsed molten rings. This value is also confirmed (dashed-dotted line) by the the relaxation 
time Rg/Dj\f for rings (diamonds). The relaxation times of linear chains can be fitted by a a power 
law with exponent 2.58 as indicated by the dashed line. 



V. CONCLUSION 

In this article we have demonstrated using the bond-fluctuation model various pro- 
nounced effects of topological constraints on the static and dynamic properties of rings. 
While the size of isolated rings scale as their linear counterparts, molten rings are quite 
compact objects characterized by an exponent u = 0.39 ± 0.03 for the chain size R. This 



is consistent with a crude Flory-like argument suggesting that one degree of freedom is 
lost for every non-concatenation constraint. 

Experimentally, melt of rings seem to show dynamical behavior very similar to linear 
chains, at least for the masses of rings usually studied |0. The same seems to be true in 
our computational study, which is of course limited to relatively modest masses - though 
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ones for which, in hnear chains, the onset of entanglements would clearly be detectable. We 
see no sign of a similar effect for rings, suggesting a larger effective entanglement length 
in the rings case. Indeed, in our simulation (which ignores hydrodynamics but respects 
the topological constraints) conformational relaxation of the rings is well described by a 
dynamical scaling involving a single characteristic time scale Tee for each chain, whereas an 
entanglement crossover would introduce a second timescale. Nonetheless, the motion is not 
that of a simple Rouse model (chains moving independently subject to local friction and 
uncorrelated noise): interchain forces are manifest in the reduced effective exponent for the 
center of mass motion at short time scales, or equivalently in an increase in the exponent 
for the dependence of the single-chain relaxation time Tee on chain length. This increase, 
fortuitously, approximately restores the exponent to the value 2 predicted by a simple Rouse 
model for gaussian rings, despite the fact that our rings are partially collapsed. 

The breakdown of a true Rouselike dynamics (which would give (73 (t) ~ t) for molten 
chains is not completely unexpected, since the Rouse model describes only one chain or ring 
in isolation. What is remarkable from the simulation is that, for the range of masses (of 
both rings and linear chains) studied so far by BFM methods, the slowing of the center of 
mass motion is consistent with a picture in which each chain or ring has to drag with it 
the neighbors in its correlation hole against the resistance of local frictional forces {D{N) ~ 
Rg{N)~^). This very simple picture may have something to offer for understanding the 
diffusion of molten linear polymers in what is classically viewed as the crossover region 
between entangled and unentangled motion, although further work is required to check 
whether the same relation holds over a range of volume fractions. Note that the diffusion 
constants and relaxation times for linear chains obtained with the BFM agree perfectly with 
simulations obtained with MD by Kremer and Grest |19| and experimental values of Richter 
by RNA PT|; these were always analyzed previously in terms of a crossover, although such 
data do not show true Rouse behavior even for small (and nor does the experimental data 
unless corrections are made for the dependence of the effective segmental mobility on chain 
length |13[). 

Likewise for rings one has to choose whether to assume a crossover or fit an intermediate 
power law to the data. The crossover interpretation is rather forced for rings, since for 
the ring sizes studied here, the intermediate (apparently power-law) behavior shows no sign 
of breaking down at either the small or large mass end of the range. Since the nature 
of entanglements in ring systems remains to be clarified, it is possible that this behavior 
could extend to quite high masses (or even, in principle, be the true limiting result for 
high A^). Certainly if an entanglement crossover is present for rings we can say that it 
is at substantially higher masses than for linear chains of the same type. In the range 
of studied, there is in particular no evidence for a crossover to a regime in which the 
rings have exponentially long relaxation times. Such behavior would anyway be surprising 
since even rings in a fixed network show an algebraic dependence of the diffusion constant 
on mass In molten linear chains, it is known that reptation cannot be the only 

mode of relaxation (the prefactors predicted from the reptation picture for the diffusion 
coefficient and for the viscosity systematically underestimate the relaxation); so there are 
probably enough alternative modes of motion to allow relaxation of molten rings on an 
algebraic timescale even in the limit of high masses. A full investigation of that limit must 
of course await further increases of computer power, especially since, as emphasized above, 



17 



any crossover to the entangled regime occurs for substantially higher masses than in the case 
of linear chains. 
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TABLE I. Mean-square bond length (6^), radius of gyration ^^^y, ring diameter (i?g) and 
diffusion constant Dn for single unknotted rings. 



N 








A 


UN 




DivlO^ 




16 


6.904 


12.9 


42.3 




12.36 


10.88 


5.2 


3.67(25) 


0.58 


32 


6.913 


25.7 


80.6 




23.30 


13.60 


6.6 


1.62(21) 


2.2 


64 


6.920 


49.3 


150.3 




42.46 


17.01 


7.3 


0.66(12) 


8.7 


128 


6.924 


92.2 


274.8 




76.66 


21.15 


7.5 


0.285(45) 


35 


256 


6.926 


169.7 


497.2 




135.36 


26.18 


6.6 


0.126(21) 


160 


512 


6.927 


304.0 


878.0 




229.88 


31.50 









TABLE II. Mean-square bond length (b'^}, radius of gyration {^Rgj, ring diameter (i?e)) ^i^S 
surface A, number of neighbors touching a reference chain n^r, monomer mobility W, diffusion 
constant Dn and rotational relaxation time Tee for unknotted, unconcatenated rings in the melt 
at volume fraction $ = 0.5. 
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